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ABSTRACT 

We use 2167 Sloan Digital Sky Survey (SDSS) quasar spectra to search for low-density oxygen in the 
Intergalactic Medium. Oxygen absorption is detected on a pixel-by-pixel basis by its correlation with 
Lya forest absorption. We have developed a novel Locally Calibrated Pixel (LCP) search method 
that uses adjacent regions of the spectrum to calibrate interlopers and spectral artifacts, which would 
otherwise limit the measurement of O VI absorption. Despite the challenges presented by searching 
for weak O VI within the Lyman-a forest in spectra of moderate resolution and signal-to-noise, we 
find a highly significant detection of absorption by oxygen at 2.7 < z < 3.2 (the null hypothesis has 
a = 80 for 9 data points). 

We interpret our results using synthetic spectra generated from a lognormal density field assuming 
a mixed quasar-galaxy photoionizing background (Haardt & Madau 2001) and that it dominates the 
ionization fraction of detected O VI. The LCP search data can be fit by a constant metallicity model 
with [0/H] = — 2.15^g Qg, but also by models in which low-density regions are unenrich ed and higher 
densit y regions have a higher metallicity. The density-dependent enrichment model by lAguirre et al.l 
(|2008f l is also an acceptable fit. All our successful models have similar mass-weighted oxygen abun- 
dance, corresponding to [(0/H)mw] = — 2.45±0.06. This result can be used to find the cosmic oxygen 
density in the Lya forest, iloxyjCM — 1.4(±0.2) x 10^^ w 3 x 10~^i7b. This is the tightest constraint 
on the mass-weighted mean oxygen abundance and the cosmic oxygen density in the Lya forest to 
date and ind icates that it conta ins ~ 16% of the total expected metal production by star formation 
up to z = 3 (jBouche et al.ll2007t) . 

Subject headings: galaxies: formation — intergalactic medium — quasars: absorption lines 



1. INTRODUCTION 

The Intergalactic Medium (IGM) is the spatially undu- 
lating reservoir of baryons from which all galaxies form. 
During the dark ages, before formation of the first stars, 
this reservoir was composed of elements in proportions 
arising from Big Bang Nucleosy thesis: hydrogen, helium 
and trace amounts of lithium and beryllium. Metals have 
been produced since the formation of the first stars and 
some of these metals were carried into the IGM. They are 
observed in the Lyman-a forest seen in quasar (QS O) ab- 
sorpti o n spectra at redshifts z < 6 (e.g. M eyer fc YorkI 
dlQSl: iCowie et al I (Il995h: iSongaila fc Cow ic (199% 
lEllison et al.l (|2000r ): ISchave et al.l (|2003f l: iPicri at jj] 
(I2006D 1. Progress has also been made reproducing 
this enrichment in models of galaxy formation with 
extragalactic winds (e.g . (Mac Low & Ferrara (1999); 
iMadau. Ferrara . fc Rees ('2001); OoDcnhcimcr & Dave 
iPieri. Martel. fc Grenon (2007); Fieri fc Martcl 

Much has been learned about the degree, extent and 
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evolution of metal enrichment of the IGM, but questions 
still remain unanswered. What is the precise level of 
enrichment in the forest, and is it consistent with that 
seen at the current epoch? Are underdense systems en- 
riched? Are all enriched regions spatially associated with 
galaxies? In this paper, we attempt to place tighter lim- 
its on the level of enrichment in the IGM as a whole, 
more specifically the low-density regime (1 ^ p/p ^ 10), 
using QSO spectra frorn the Sloan Digital Sky Survey 
(SDSS. lYirkeiTall 120001) . We have searched for O Via 
(A1032A) at redshifts 2.7 < z < 3.2. At the densities 
analyzed, O Via is expected to be the strongest heavy el- 
ement transition in the optical band ([Ranch et al.lll99"7l : 
iHellsten et al.l [r99^. However, the line happens to re- 
side in the same wavelength range as the Lyman series 
lines at other redshifts resulting in severe confusion as 
to the identity of the absorption. Any search for O VI 
must involve a sophisticated approach to overcoming this 
problem. 

Fitting of narrow lines or lines that also show 
Lya a bsorption have been performed by various au- 
thors (ICarswell et a l.' '20021 jsimcoe et alj l2002l 12004 
[Bergeron fc Herbert- Fort 2003) order to find individ- 
ual O VI systems, where high resolution spectra are 
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rected toward finding metals distributed throughout the 
forest without identifying individual systems. In the case 
of the search for O VI, they are limited by cosmic vari- 
ance in the contaminating Lyman lines and noise. The 
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standard pixel (SP) correlation search (so metimes known 
as 'p i xel optical depth ' ) for O VI (e.g. lApruirre et alJ 
I2002t iPieri fc Haehneltl [20041: lAguirre et al.ll2008l here- 
after AOS) is not, however, a suitable approach for the 
detection of metals in the SDSS Lya forest because of the 
moderate resolution {R — A/AA = 1800) and signal-to- 
noise ratio. A pixel correlation search is possible in prin- 
ciple, but modifications to the approach are required. We 
have developed a tailored pixel method, which we refer 
to as the Locally Calibrated Pixel (LCP) search method. 
Our statistical analysis of O VI absorption in the low den- 
sity IGM complements the direc t search for strong O VI 
absorption in the same spectra ([Frank et al.ll2009al lbl). 

While interpretation of pixel searches using high res- 
olution, high S/N spectra is more straightforward, the 
small amount of high quality data leaves the conclusions 
subject to uncertainty owing to cosmic variance. The use 
of the large number of SDSS spectra (> 2000 used here) 
with the LCP search method offers a large reduction in 
this variance (and potentially an increase in overall S/N). 

All pixel correlation searches require simulations in or- 
der to ascertain the gas mass densities of enriched sys- 
tems, the number densities of ions required, the ioniza- 
tion fraction (and so metallicity of those systems), and 
the impact of observational limitations and errors. This 
is particularly important here as the resolution of the 
SDSS spectra introduces a significant uncertainty in the 
gas mass density of enriched systems. To address this 
issue and interpret our search results, we have simulated 
each SDSS spectrum used 10 times for various models of 
oxygen enrichment, totaling more than 4.5 million syn- 
thetic SDSS spectra. 

Noise levels in individual SDSS spectra are higher than 
those of high resolution spectra used in previous studies, 
and with this (along with the lower resolution) comes 
more uncertainty in the continuum fitting. This prob- 
lem becomes increasingly severe at the blue end of the 
spectrograph, where most of the signal of O VI lies. The 
correlation search technique that we outline below is cal- 
ibrated to the flux level in the local spectrum to take into 
account variations in the continuum (hence the name 'lo- 
cally calibrated pixel search'). 

This paper is set out as follows: in ^j2|we describe the 
data set used, in ^J3|we describe the LCP search method 
and results, in fJHwe describe our interpretation of these 
results using simulated spectra, and in Sj5]we describe our 
measure of mass-weighted mean oxygen abundance and 
the cosmic oxygen density. This is followed by a discus- 
sion of the ionization corrections and other searches and, 
finally, the conclusions. 

Following standard notation, we write 
[0/H] =log(0/H)-log(0/H)o, where 0/H is the 
ratio of oxygen to hydrogen (by number). We adopt 
log(0/H)(7) = —3.13 e stimated for the solar envelope by 
iGrevesse et al.l ([1996( ). for ease of comparison with other 
studies of oxygen enrichment of the IGM. While some 
analyses based on more detailed ste llar atmosphere mod- 
els have argued fo r lower (O/H)© ([Asplund et al.|[2005[ : 
[Caffau et al.|[2008l) , stellar interior models and helioseis- 
mology strongly supp o rt the Grevesse et al. abundance 
scale ([Basu fc Antia[ [20041 : iDelahave fc Pinsonneault[ 
20061). 



For values of solar abundance taken from 



Asplund et all ([2005D and [Caffau et al.l ([2008[ ). scale our 



[0/H] results by -1-0.21 and -1-0.11 respectively. 



2. THE SAMPLE OF SPECTRA 

The O VI doublet starts to become visible at the blue 
end of SDSS quasar spectra at redshifts of ^abs > 2.7. 
Hence, we retrieved QSO spectra from the SDSS Third 
Data Release (DR3; [Abazaiian et al.ir2005[ ) beyond this 
lower boundary redshift using the QS O Absorption Lin e 
Systems (QSOALS) project database ([York et al.[[2006[ ). 

The spectra we used were fitted with an automatic con- 
tinuum procedure developed by Arlin Crotts and used in 
the QSO absorpt ion line data base of the SDSS QSOs 
([York et al.|[2006[ ) . The continuum is derived by discard- 
ing absorption-contaminated pixels until the variance in 
the remaining pixels corresponds to that prescribed by 
the noise model for the measurement. This provides a 
smoothed estimate of the spectrum, forming a "contin- 
uum" for the measurement of narrow lines over a range 
of pixels that varies with the large scale variations in the 
SDSS spectrum. At a minimum (e.g. near narrow emis- 
sion lies from the QSO or the sky, or at the steep drop at 
the edge of a broad absorption line) the smoothing is over 
a 20 pixel region (between IsA to 24A in the forest used 
here). In smooth regions of the spectrum, outside the 
Lyman-a forest, the method produces continua accurate 
to about 2% when the photon noise is small enough to 
permit such precision. The method produces a reason- 
able estimate for the continuum in the Lyman-a forest 
for our purposes. Note that the rare cases of a Lyman-a 
BAL or the more common cases of an O VI BAL will be 
fitted away and this is a desirable feature for our pur- 
poses. 

There are rarely rapid continuum changes to be ex- 
pected over 30 A, and certainly no significant non-linear 
trends. This gives us a good enough continuum for ap- 
plication of a locally calibrated continuum for the O VI 
absorption measurement, as described in The contin- 
uum fitting method may subtract large-scale structures 
in Lya absorption on scales > 2800 kms""'^, but this ef- 
fect will be minimal as indicated by the good agreement 
between the probability distribution function (PDF) of 
the Lya optical depth in the data and the models as 
shown in Figure |4| below. This continuum fitting method 
results in a ^ 16% systematic underestimate of the con- 
tinuum (and so an overestimate of the transmitted flux) , 
measured as an equivalent underestimate of the mean 
flux decrement in the sample with respect to values in th e 
literature (e.g [McDonald et al.[ [20061: [Kim et al.[ [20071) . 
We have taken this into account by introducing the same 
systematic to our simulations as described in ^4.11 

We do not exclude BAL quasars from the sample used 
in this study. We have tested the impact of these BAL 
QSOs by performing the L CP search excluding the 430 
QSOs that are identified bv [Gibson et al.l (|2009[ i as BAL 
QSOs and fall into our sample. The impact of this mod- 
ification is much smaller that t he la error bars. 

As can be seen in Figure 3 of [Frank et all ([2009a[ ). the 
majority of the sources have i band magnitudes between 
19.5 and 20.5. The average signal-to-noise ratio of the 
spectra within the region of interest for us is only about 
2.5, much less than the overall signal-to-noise ratio of 
the spectrum, since the flux levels within the Lyman-a 
forest are much lower than those beyond it. This fig- 
ure also shows the distribution of these average signal- 
to-noise ratio values with the emission redshifts of the 
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Noise characteristics at 5000 
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Fig. 1. — Example of the noise characteristics of the sample 
spectra at 5000 A. Plotted in black are the estimated noise levels 
in real flux, (j{f\), versus the flux at 5000 A. The noise can be 
separated into two components, the detector dependent read-out 
noise, Uro, and the flux dependent photon noise, rj^^. The dotted 
red line shows the two parameter flt to the data, as described in 
the text. The dashed black lines indicate the location of the noise 
if it were purely dominated by the Poisson statistics of the photon 
noise with 5, 10 and 20 counts. Where the photon noise dominates, 
the counts are clearly above 20, and so this Poisson distribution 
is well described by a Gaussian appoximation. Elence the noise 
distribution is Gaussian at all levels of flux. We have confirmed 
this for the complete wavelength range covered by SDSS. 

sample sources. Particularly towards higher redshifts, it 
is apparent that the increasing density of the forest and 
the general faintness of the sources lead to a severe de- 
crease in the signal-to-noise ratio. 

In order to ensure that the noise characteristics of 
the sample can be adequately reproduced in the simu- 
lated spectra, we have performed an analysis of the be- 
haviour of the noise in the full data set at every observed 
wavelength. The uncertainty in the flux measurement, 
o'I/a, A), is the sum in quadrature of two terms, the read- 
out noise of the detector system and the photon noise. 
The former is a Gaussian distribution with (7^0 (A) = fci, 
where fci is a constant at a particular wavelength, in- 
dependent of the incident flux. The latter is a Poisson 
distribution with crph{f\,X) — k2 x vTa, where k2 is a 
constant at a particular wavelength. 

At each wavelength, we fit the distribution of the error 
estimate in the flux provided by the QSOALS database 
with this two parameter model, and hence are able to 
describe the complete noise characteristics of the under- 
lying SDSS sample by the wavelength dependent A:i(A) 
and fc2(A). Figure [T] shows an example of the noise dis- 
tribution as a function of the estimated flux, fx, at a 
specific wavelength. We have used this formalization of 
the noise characteristics to produce a realistic represen- 
tation of observing conditions in synthetic spectra as de- 
scribed in i j4.1l This figure also illustrates how the noise 
may always be characterized as a Gaussian distribution, 
since the (Gaussian) read out noise dominates where the 
Poisson nature of the photon noise is significant. 

3. THE LOCALLY CALIBRATED PIXEL SEARCH 

Pixel correlation searches for the detection and inves- 
tigation of metal enrichment in the IGM have, until now, 
been limited to high resolution quasar absorption spec- 



tra. In these spectra, the optical depth of every pixel of 
the Lyman-a forest can be characterized as a measure 
of the density in the IGM. Where the Lya absorption is 
saturated (and so the measured Lya optical depth mea- 
sure is noisy), higher order Lyman lines can be used to 
reconstruct the Lya optical depth. 

Ionization species other than neutral hydrogen can pro- 
duce their own forest of absorption. In particular there is 
a forest of five-times-ionized oxygen (O VI) and a forest 
of three-times-ionized carbon (C IV). These metal line 
forests can also be dealt with on a pixel-by-pixel basis 
at each redshift. Hence, one can collect Lya and O VI 
or C IV pixel pairs, bin by Lya optical depth, and one 
would see a clear trend in metal line optical depth (us- 
ing only basic assumptions of ionization and metallicity, 
which we will return to throughout this paper). 

Difficulties arise when the measured optical depth of 
the metal pixels is not entirely due to the intended 
species. Various techniques have been used in the SP 
search to minimize the degree of contaminating absorp- 
tion, but it cannot be removed entirely, and it tends to 
dominate where absorption is weakest. These techniques 
involve correction of contaminating Lyman series lines 
in the O VI region of the spectra by removing expected 
higher order lines (based on the Lya absorption at the 
same redshift and longer wavelength), or the use of both 
O VI lines to derive a minimally contaminated signal. 

As is done in all pixel correlation searches in quasar 
spectra, we compared the absorption from a species at a 
fixed redshift with the absorption from another species 
at the same redshift. Where the absorption is corre- 
lated, a detection is found. In this analysis we search 
for O Via (AIO32A) absorption correlated with Lyman- 
a (A1216A) absorption and make no use of the O VI 
doublet line (O VIb) at A1038A. We produce these pixel 
pairs for a range of redshifts and on a grid set by the 
spectral binning of the O Via pixels (the Lya absorp- 
tion is interpolated). Nearby pixels with a wavelength 
close to the O Via wavelength are used to characterize 
continuum errors and contaminating absorption. This 
procedure is explained in more detail below. 

We have simplified the algorithm in recognition of the 
fact that there are features of the SP search that are un- 
suitable in the context of our local calibration method. 
These techniques are directed toward subtracting con- 
tamination. We do not seek to remove contaminating 
absorption; as described above, we characterize it, and, 
as shall be seen below, the LCP analysis does this auto- 
matically. 

The lower redshift limit for every quasar in the sample 
is set by the blue end of the spectrograph (around 3800A) 
or the Lyman-7 line at the quasar redshift, whichever re- 
sults in a higher minimum redshift. We discard pixels 
with Lyman-7 as they provide little extra information 
but modify the distribution of absorbers in O VI sig- 
nal enough to introduce unnecessary uncertainty. Note 
that this requirement also means that no absorption 
at the Lyman limit is possible. We rule out regions 
within 5000 kms^^ of the QSO emission redshift to elim- 
inate m ost effects due to the QSO and its environment 
(|Wild e t aL 200S). This sets the upper redshift limit and 
also results in an O VI signal entirely within the Lyman- 
[3 forest (since the separation between O Via and Ly/? 
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Fig. 2. — Redshift distribution of pixel pairs available for the 
analysis in the SDSS sample. The lower envelope is set by the 
lowest redshift systems observable in O VI at the blue end of the 
spectrograph. The upper envelope is set by 5000 kms~^ below the 
emission redshift of QSOs in the sample 

absorption at the emission redshift is 1,804 kms"^ in 
the O VI frame) . 

We combine the pixel pairs for each QSO into one com- 
bined sample for all the SDSS spectra in the required 
redshift range producing over a million pixel pairs. The 
redshift distribution of these data is shown in Figured 

3.1. H I Absorption 

The determination of the degree of H I absorption is 
largely unchanged from the SF search. The optical depth 
to Lyman-a absorption Ti^ya is determined for each pixel. 
We make the same noise requirement as other searches: 
the flux transmission must be crnoiso/2 from either the 
zero level or the continuum, where (Tnoise is the contin- 
uum scaled error estimation in the continuum fitted spec- 
trum. This is a rather strict criterion in the context of 
the SDSS sample and discards « 20% of pixel pairs. We 
do not attempt to extend the analysis to the higher den- 
sity regime, where Lya absorption is saturated, by using 
higher order Lyman lines to derive the Lya optical depth. 
This is a regime to w hich other methods are better suited 
([Frank et al.ll2009allbl l . It should also be noted we assume 
that all O VI is associated with H I absorption, but under 
most physical conditions this is a fair assumption. 

The large amount of data at our disposal allows us to 
aggressively discard less reliable measures of TLya • Hence 
we require that Cnoise < 0.1, which discards a further 
70% of the pixel pairs, so the feature we are searching 
by (Lya) is as clean a signal as we can make it while 
preserving a large sample of data. As we shall show, 
the feature we are searching for (O Via) is treated with 
no such noise requirements, and brute force averaging is 
used to provide a clean signal. 

3.2. Measuring the Bulk O VI Absorption 

Any measure of weak correlated O VI absorption must 
deal with the biggest uncertainty: the continuum level. 
The primary goal for the following technique is to re- 
move the impact of errors in the continuum fitting from 
the derived signal. Rather than simply taking the trans- 
mitted flux at the position in the spectra corresponding 
to the required O VI absorption we determine the ratio 
of the transmitted flux at this point compared to that 
of a neighboring pixel. Where this pixel is sufficiently 



close, this gives us a measure of the absorption with the 
local error in the continuum fitting subtracted. It also 
results in a subtraction of average levels of contaminating 
Lya, Ly/? and O VIb absorption along with a measure of 
nearby O Via absorption. The ratio of transmitted flux, 
F, at the O Via pixel to that on a locally calibrating pixel 
is 

-Fovia _ fc exp[-rovia] exp[-T^] -I- noise 

Fi fed exp[-Tovia,;] exp[-Tu,i] + noise; ' 

where fc is the flux level at the continuum, Tovia is the 
optical depth to O Via absorption, 'noise' is the Gaus- 
sian instrumental noise, the subscript T indicates a local 
value and the subscript 'u' refers to uncorrelated con- 
taminating absorption. This uncorrelated absorption is 
dominated by Lya and Ly/3 absorption, but it also in- 
cludes contamination by the second member of the O VI 
doublet (tu = TLya -|- TLyi3 + Tovib), all at different red- 
shifts, roviaj is the local O Via absorption, which is also 
uncorrelated where the appropriate local pixels are used, 
as we shall argue in ^'3A\ 

Multiple local pixels are used, and each one provides 
a new measure of spectral characteristics local to our 
required O Via pixel. Operationally, we treat each local 
pixel as providing a new measure of i^oviaZ-F) and so a 
new pixel pair of TLya — ^bvia/^/- In this way, the final 
statistical error is dominated by the number of Lya - 
O VI pixels and not by the noise in the local pixels. 

3.3. The Aggregation of Pixel Pairs 

Since we have a large sample of spectra with compara- 
ble characteristics and a limited redshift range we aggre- 
gate all acceptable pixel pairs from all our spectra. We 
use pixels in a redshift range of 2.7 < z < 3.2 both to 
keep our data as homogeneous as possible and to limit 
ourselves to the subset that provides the strongest signal. 
This limits us to 2167 QSO spectra in our sample, out of 
the 5767 z > 2.3 quasar spectra in SDSS DR3. 

Once we have a list of TLya — Fovia/Fi pairs, as set out 
above, we bin them by their TLya and take the median of 
the TLya and the -Fbvia/f; in order to produce our final 
search for O Via. We use the median as an outlier- and 
noise-resistant measure of the typical absorption. 

Since r„ ^ and t„ are drawn randomly from the same 
distribution, and (where the local pixel is sufficiently 
close) fed ~ fc, the median(_Fovia/-F)) is a quantity that 
characterizes correlated O Via absorption normalized to 
a factor dependent on the mean O Via optical depth. 
The precise level of this normalization is unimportant 
in the following results; as with other pixel correlation 
searches, a statistically significant correlation is required 
for a detection of metals. It is, however, broadly indica- 
tive of the average O Via absorption level. It is also 
notable that this level is recovered in the null tests that 
follow. 

3.4. Choice of Local Pixels 

For local pixels we chose those that are i) sufficiently 
far away from the O Via search pixel to avoid comparison 
with pixels within the same extended complex associated 
with the same enriched region, ii) sufficiently close to 
avoid producing spurious correlations with the Lyman-/? 
and O VIb absorption, and iii) sufficiently close to mea- 
sure the same continuum level. The second requirement 
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ensures that the third requirement is satisfied. Lyman- 
l3 is stronger and correlated over a larger velocity range 
compared with O VIb, hence we design our choice of lo- 
cal pixels for avoidance of this line. We also require that 
the choice of local pixels are mirrored on each side of 
the search pixel in order to cancel continuum level error 
trends across the line. 

We take 18 local pixels in total for this analysis: the 
9th to the 17th pixel from the search pixel at both lower 
and higher wavelengths. This results in separations of 
> 600 kms~i and < 1200 kms^^ from the O Via and 
Ly/3 transitions, which is approximately 9 — 16A in the 
observed frame at z = 3. This provides « 800, 000 pixel 
pairs in all the searches that follow. 

We tested our choices of locally calibrating pixels by 
treating each of them individually as a 'search target' 
(searching for correlations with Lya) and finding the ex- 
pected null result. Pixels within 600 kms~^ of O Via or 
Ly/3, on the other hand, yielded signals of contaminat- 
ing correlated absorption. The 600 kms^^ scale is char- 
acteristic of the winds from Lyman break galaxies (e.g. 
iPettini et aLll2001l:lAdelberger et al.ll2003D and the max- 
imum clust e ring scale found for C IV (lAdelberger et al.l 
[200l [20051 iPieri at al.l [200l iScannapieco et al.l 120061 )7 
so these results seem reasonable for O VI. Correla- 
tions in Ly/3 extend o ver greater scales than this due 
to large scale structure (jMcDonald et al.ll206 6). but such 
lines are somewhat stronger than expected here and our 
test searches confirm this. In principle C II(A1036A) 
could produce correlated absorption within our calibrat- 
ing pixel range, but in practice it appears to be too weak. 

In test searches we detected the doublet line O VIb 
at a consistent level with that of the O Via line. Since 
this line merely constitutes a parallel but less sensitive 
measure of O VI that does not contribute to the search 
algorithm, we have not used this information. 

3.5. VI Search Results 

Once the pixel pairs required are combined and binned 
by their TLya, we find their medians to plot Foyi/Fi vs. 
Thya as shown in Figure |31 We required that more than 
200 pixel pairs be available for a bin to be used. In all 
panels the error bars are produced by bootstrap resam- 
pling with 100 realizations and a bootstrap element of 1 
SDSS QSO. The median is found for every realization, 
and the error bars are then taken to be the la deviation 
in the distribution of medians. 

Figure El shows Fovi/Fi vs. TLya using pixels in the 
redshift range 2.7 < z < 3.2. There is a clear correlation 
between TLya and Fqw/Fi, and this correlation extends 
to all measured values of TLya . There is both excess O VI 
absorption in pixels with strong Lya absorption and a 
deficit of O VI absorption {Fovi/Fi > 1) in pixels where 
Lya absorption is weak. 

The null test shown in Figure [3] is performed using 
the same technique as set out above with the only mod- 
ification that the restframe wavelength for correlations 
is no longer the O Via restframe wavelength (A1032A) 
but A = 1047A. As can be seen, the search approach 
passes the null test by recovering the expected result of 
^1047 ~ Fi, where neither an excess nor a deficit in ab- 
sorption is seen as no signal is present. 

Figure |3] is the main observational result of this pa- 
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Fig. 3.— The LCP search for O Via in the SDSS Lya forest 
from pixels with 2.7 < z < 3.2. The red filled circles are the median 
F/ Fi on and off the 1032A O Via resonance line vs the TLy^ at the 
same redshift. Where Fqyi/Fi is below (above) 1 there is a greater 
(lower) than average degree of O VI absorption. A correlation in 
this plot is indicative of a detection of O VI. The black '-'s are the 
null case of searching on and off 1047A. 

per. In high resolution spectra there is an approximate 
one- to-one relation between tlv q and gas o verdensity 
o/o llRauch et al.lll997l : [Croft et a l. 1998; Wein berg et"al] 
11998 ). At SDSS resolution this is no longer the case, so 
the interpretation of this result in terms of enrichment of 
the IGM must rely on synthetic spectra drawn from an 
underlying physical model. 

4. INTERPRETATION USING SYNTHETIC SPECTRA 

4.1. Production of Synthetic Spectra 

We use an approach for the production of simulated 
line-of-sight dens i ty dis tributions that was set out in 
IPieri fc Haehneltl (|2004D . We restate this method here, 
but for a full description refer to that paper. We 
use co smological parameters consistent with WMAP5 
(jKomatsu et al. 2008) {Qm = 0.268, Qa = 0.732, = 
0.0441,(78 = 0.776 and Ho = 70.4 kms-^Mpc^i) in these 
simulations 

We start with a 3D power s pectrum of dark matte r 
density fluctuations taken from l Efstathiou et al.l (|1992D . 
This is normalized using the cosmological parameter, tJg 
(the amplitude of clustering on an 8h~^Mpc scale at the 
present epoch) , and the linear growth of structure at the 
required redshift. A Jeans length filtering (~ 0.9ft.~^Mpc 
comoving in the surveyed redshift range) to remove struc- 
ture on small scales is used to describe baryonic pressure 
effects. 

Two ID power spectra are produced by integrat ion of 
the 3D power spectrum (jKaiser fc Peacocfl |1991D . and 
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Gaussian random realizations of structure are generated. 
These ID realizations of structure are combined to pro- 
duce coupled density contrast and velocity fields. This 
velocity field, v{x), is a non-dynamical approximation 
based on the linear de nsity field, S{x). For m ore details 
on this approach see iBi fc Da vidsenl ()1997( ) and refer- 
ences therein. 

Further following the method of iBi fc Davidsenl (|1997l ). 
we convert our linear density field to a lognormal distri- 
bution, using the mapping 



p/p{x)=eMS{x)~{S^)/2], 



(2) 



which they find reproduces the probability distribution 
function of non-linear structures. It should be noted that 
this method fails to reproduce the clustering of non-linear 
structures (Viel et al. 2002). The density field produced 
has a resolution of O.OlMpc comoving, which substan- 
tially over-samples the resolution required by the SDSS 
sample but is adopted in order to ensure that the physics 
of the IGM is well described. 

The temperature of the IGM in these simulations is 
set by the balance between photoionization heating and 
adiabatic cooling and follows the power-law relation 

T^Toip/p)'-^ (3) 

(|Hui fc Gnedin|[l997l ). where p/p is the overdensity and 
the temperature at mean density is Tq = 2 x 10'*K. The 
temperature is required to calculate the Doppler param- 
eter for this gas. 

For a given set of cosmological parameters and a he- 
lium fraction of 0.242 by mass, it is trivial to calcu- 
late the hydrogen density. Oxygen is added in quanti- 
ties as described in the following section. The ioniza- 
tion fractions of hydrogen ai id oxygen are calcu lated us- 
ing CLOUDY ver s ion 08 .00 (|Ferland et allll998H and the 
iHaardt fc Madaul (|2001f ) (hereafter HMOl) quasar and 
galaxy UV background model (designated QG) with a 
10% escape fraction from galaxies. While these calcula- 
tions do take into account collisional ionization, the lack 
of heating other than via photoionization in our simula- 
tions makes the conditions for it to occur rare. As such, 
we predominately interpret the detectable O VI to be 
outside of this coUisionally ionized, hot gas phase. This 
assumption is supported by observations that the bulk 
of previously detected O VI lines in the Lyg forest are 
photoionized (ICarswell et al.ll2002t [Bergeron et al.ll2002l: 
ISimcoe et al.ll200 4t) based on their line widths. This ap- 
proach provides fully self-consistent models since regions 
with p/p < 10 will be overwhelmingly photoionized and 
we find this regime dominates our measured O VI (see 

We calculate the optical depths for Lya, Ly/3, O Via 
and O VIb using 



t{v) 



Tre^A/ f n{v') 



eCH. 



b{v') 



exp 



\ b{v') 



dv', (4) 



where the transition dependent parameters are the rest 
frame wavelength A, the oscillator strength /, the num- 
ber density of absorbing particles n, and the Doppler 
parameter b. Hz is the redshift dependent Hubble pa- 
rameter, V is the line-of-sight velocity, Vq is the velocity 
of the Hubble fiow, e is the electron charge, nif. is the 
electron mass, and c is the speed of light. This is inte- 
grated over the narrow range of velocities that contribute 
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Fig. 4. — The PDF of Lya optical depths in the observed sample 
{solid line) and one simulated suite {dashed line). The simulated 
sample was produced with a quasar and galaxy HMOl UV back- 
ground, the intensity of which was renormalized to match the ob- 
served distribution giving an ionization rate of F = 13 X lO^^'^s"^. 

significantly to the calculation (effectively a few times the 
Doppler parameter). 

We add instrumental broadening to match the SDSS 
spectral resolution. As stated in Sj2l we have performed 
an analysis of the characteristic a{fx,X) in our sample, 
and we use this to add Gaussian noise to the spectra. In 
order to achieve this we use the continuum level in the 
real spectra to rescale the simulated transmission fiux to 
a fiux, f\. The photon noise is well described by a Gaus- 
sian distribution where it dominates, and the readout 
noise is also a Gaussian, so treating the combined noise 
f (/a, A) as a Gaussian is a good approximation Fi- 
nally, we rescale our transmitted fiuxes up by 16% in line 
with the systematic underestimate of the continuum in 
our sample 

The intensity of the UV background was renormalized 
(see H6.ip to provide a match of the PDF of optical depths 
derived from our synthetic spectra with that of the ob- 
served sample. We obtain agreement with an ionization 
rate of F = 13 x 10~^^s~^ as is shown in Figure |4l The 
quality of this agreement is sufficiently good to populate 
the bins in TLya with statistics that are representative of 
the LCP search in the data. 

4.2. The Monte Carlo Approach 

For every oxygen abundance pattern we simulate a full 
suite of 2167 SDSS QSO spectra and do this ten times. 
For each one we perform an LCP search for O VI, and 
the resultant Fovi/Fi vs. TLya is averaged over the ten 
suites to obtain a model prediction (so that the quality of 
the fit between the model and observed LCP searches is 
dominated by the error in the observed search) . We have 
used this set of 21,670 line-of-sight density distributions 
with varying levels of oxygen abundance to interpret our 
observations. Each of these lines-of-sight require the sim- 
ulation of a ID density distribution of length around 
lOOO/i^^Mpc. The production of this quantity of data 
is not practical with fully numerical simulation methods 
given current computer processor limitations. This un- 
derlines the value of the analytical technique described 
in the previous section. 212 models of oxygen enrich- 
ment have been tested, totaling more than 4.5 million 
synthetic spectra produced. 

In the following analysis we do not include the mea- 
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Fig. 5. — The LCP search compared with models of density 
independent [O/H]. The red filled circles are the LCP search as in 
Figure \3\ Four models are shown: the solid line is the null case 
where no oxygen is added, the dashed line shows the best fitting 
case of a constant [0/H]= —2.15. A further two models show 
oxygen metallicities that are ruled out: [0/H]= —2.00 {dot-dashed 
line) and [0/H]= -2.30 (dotted line). Note that the TLy^ ~ 3.3 
point is not included in the fits. We draw detailed conclusions on 
metallicity constraints in JS] 

sure at T^ya ~ 3.3 (see Figure |3|). The inclusion of this 
point results in poor agreement with the models. This 
is consistent with the findings of AOS, that high TLya 
absorption results from a different population of O VI 
absorbers that are most likely coUisionally ionized. 

4.3. Comparison of Real and Synthetic Spectra 

We add oxygen with a constant [0/H] to our synthetic 
spectra and have varied this [0/H] to find a good fit to 
the observed search. Figure [5] shows a comparison of the 
observed search (filled circles) with four models. Once 
more we see that our LCP search responds well to the 
null test (solid line) by returning a result of Fqvi ~ Fi 
when no oxygen is added (the reduced — 10). In 
the best fitting case, [0/H] = —2.15 (dashed line), we 
find good agreement with the observed result, and there 
is no need for the addition of more parameters in the 
fit. Models with constant [0/H] = -2.0 (dotted line) 
or [0/H] = —2.3 (dot-dashed line) are clearly ruled out 
and indicate that the LCP search in the SDSS sample is 
sensitive to small changes in the abundance of oxygen in 
the IGM. 

We can quantify the quality of the fit by using the 
test where the error bars are given by the boot-strapping 
of the observed data. We have done this for varying 
[0/H] as shown in Figure [6] (solid line). The 1.5(7 and 
2(7 level are shown as the dotted lines. It is clear that 
our best fitting model with [0/H] = —2.15 constitutes 
an adequate fit to the observed search. We can give this 
measure a 2(7 confidence interval of [0/H] — — 2.15j^Q Qg. 
These models constitute the range of constant [0/H] 
models favored by the LCP search. The quality of this 
fit is suggestive of an underlying density dependence to 
the [0/H] signal (as found by AOS), but the statistics do 
not warrant a two parameter fit. 

We may not be able to constrain the density depen- 
dence of this oxygen enrichment, but we can place limits 
on the density range in which it resides. Hence we in- 
troduce a new parameter: a gas overdensity threshold, 
(p/p)cut- Where plp< {plp)cat, we add no oxygen, and 
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Fig. 6. — The test of fit [solid line) between LCP searches per- 
formed on our models with varying density independent [O/H] and 
our observed LCP search. The dotted lines show the la, 1.5cr and 
2(T confidence levels indicating the range of [O/H] that constitutes 
a good fit. 

where p/p > {p/p)cut we add the constant [0/H] spec- 
ified. Figure [7] shows a likelihood plot of the fit to the 
SDSS data set for 20S models (with two parameters and 
nine data points). Where (p/p)cut = 0, no threshold is 
used and the models correspond to those used in Figure 
H There is little or no change in the quality of the fit 
going from {p/p)cut = to (p/p)cut = 1, and this indi- 
cates that our LCP search for O VI in the SDSS QSO 
sample is not sensitive to metals in underdense regions 
of the Universe. 

As the threshold is raised in our models a degen- 
eracy becomes apparent, as seen in Figure [T] remov- 
ing metals from low-density regions requires one to add 
more metals to higher density regions to compensate, 
and it is still possible to obtain a good fit. This indi- 
cates that our search is rather insensitive to the den- 
sity of the systems enriched. As expected, the tight 
relationship betwe en TT,vr» and overdensity seen in high 
resolution spectra (iRauch et all 119971: ICroft et aLlll99Sl: 
[Weinberg et al.lll99SD breaks down due to the moder- 
ate resolution of the SDSS sample. Dense systems can 
always produce a low apparent TLya by virtue of the in- 
strumental broadening over multiple pixels, but the key 
question is to what extent those systems can reproduce the 
observed signal? The agreement begins to break down at 
(p/p)cut ~ 6, where we have 13% confidence {1.5a) in 
our fit to the observed result, and it breaks down com- 
pletely at (p/p)cut ~ S.5, where we have a 5% confidence 
(2(7) in our fit. This breakdown occurs because there are 
proportionately too few higher density systems to obtain 
a fit to the search curve even with arbitrarily high metal- 
licity. The best fit is obtained with {p/ p)cut = 4, where 
the fit has — 9-4 (23% confidence for two parameters 
and nine data points). 

While Figure [7] shows that a range of two-parameter 
models can fit the data, we will show in fJSlthat the mass- 
weighted mean oxygen abundance is well constrained. To 
demonstrate this, we must first address the issue of scat- 
ter in the metallicity. 

4.4. Scatter in the Metallicity 

ISchave et al] ()2003[ ) (hereafter S03) and AOS found 
that a lognormal scatter in the metallicity is a good fit 
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Fig. 7. — The l.Scr (dark red) and 2a (orange) likelihood of fit 
between LCP searches performed on our models and our observed 
LCP search (with two parameters and nine data points). The 
models have a constant [O/H] above a threshold gas overdensity, 
(p/p)cut. Below this threshold no oxygen is added. The contours 
characterize the range of densities for which this search for O VI is 
sensitive. The likelihood is largely unaffected by the addition of a 
threshold at p/pcut = Ij indicating that the search is not sensitive 
to metals in underdense systems. Where a good fit can no longer 
be obtained for a particular value of (p/p)cut, this constitutes the 
upper limit on the density range to which we are sensitive. 

to all percentiles in their search. We assume the same 
scatter in order to have a better understanding of the 
nature of our search results. In particular, it is not clear 
if an [O /H] model without scatter that fits the data gives 
a measure of the mean or median metallicity in the IGM. 

A model with lognormal metallicity scatter, i.e Gaus- 
sian scatter with variance in \ogZ = [O/H], has 

\ogZmed^median{\ogZ) = (logZ) 

= log(Z)-i^af, (5) 

where logZ^ed and (Z) are the median and mean metal- 
licity , respectively, and we adopt cr^ = 0.5 (jSchave et al.l 
l2003l eq. 7). The lognormal scatter is added on 
scales of 1.2/i~^Mpc comoving in regions where —0.5 < 
\og{p/p) < 2. As equation [5] shows, lognormal models 
with different scatter, as, have the same (logZ) and the 
same Z^ed but different (Z). 

In order to test whether the LCP search for O VI is 
sensitive to the mean or median oxygen abundance we 
explore both interpretations. Figure [5] compares our best 
fitting constant [O/H] model with no scatter (as shown 
in Figure [5]) to two models with scatter, one that has 
the same {Z) and one that has the same Zmed- The 
good agreement in the former case indicates that our 
search for O VI is sensitive to the (Z) or the mean O/H. 
Furthermore, we have tested this hypothesis using the 
A08 density-dependent fit to O VI seen in high resolution 
spectra. They report a model fit 



[O/H] = [0/C] + [C/H] 



3) 



-0.65 



+0.1 



{\og[p/p] - 0.5) 



(6) 



that is derived from O VI correlated with C IV and C IV 
correlated with H I. This model is a three parameter fit 
to the median [O/H]. Since it is fit to independent data 
we treat it as having nine degrees of freedom here (equal 
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Fig. 8. — The impact of scatter in our models and its conse- 
quences for the nature of the LCP search. The red filled circles 
are the LCP search as in Figure \3\ The dashed line is the best 
fitting model with constant [O /H] and no scatter as in Figure [5] 
We have added lognormal scatter to this model while keeping ei- 
ther the mean O/H fixed {dotted line) or the median O/H fixed 
{dot- dashed line). The agreement of the former but not the latter 
indicates that the scatter-free model fit yields a robust estimate of 
the mean O/H, not the median O/H (or the mean [O/H]). Note 
that the Tj^ya ~ 3.3 point is not included in the fits. 

to the number of data points). When we reproduce this 
model in our simulations including scatter and perform 
the LCP search, it provides a reasonable fit to the SDSS 
data (about 2(t). However, where we have simply added 
metals to our simulation at this median level without 
scatter our agreement is poor (reduced = 6). When 
we use the correction factor in equation [5] and add metals 
at the mean level instead, we obtain a reasonable fit once 
more (at about 1.5(7 confidence). 

5. MEASURING THE MEAN OXYGEN ABUNDANCE AND 
THE COSMIC OXYGEN DENSITY 



Sections 14.31 and 14.41 present a variety of models that 
are consistent with our LCP search results, despite dif- 
ferences in density dependence and scatter. We can ask 
whether these models have a common feature of oxygen 
abundance that is well determined by our data. To this 
end we have calculated the volume-weighted and mass- 
weighted means of O/H for all models that fit the data: 
those within 1.5a and 2a contours, along with the AOS 
model. Specifically, we compute the means by summing 
O/H (not [O/H]) over all pixels in the simulated spec- 
tra with mass density in the range —0.5 < log[p/p] < 2, 
either weighting each pixel equally (volume-weighted) or 
weighting in proportion to mass density. 

The mass- weighted mean oxygen abundance, 
[(0/H)mw]j for the 'ridge line' of best fitting mod- 
els at every {p/p)cut is shown in Table [TJ We list the 
models without scatter, but models with scatter yield 
the same results since (as shown in Figure [8|) they must 
have the same mean oxygen abundance to fit the data. 
Taking models within the 1.5 and 2a contours (rather 
than just this ridge line) provides a measure of both 
uncertainty in density dependence and observational 
uncertainty. 

While the volume- weighted abundances span 0.4-1.3 
dex in [O/H], the successful models all have mass- 
weighted means in the range —2.51 < [(0/H)mw] < 
-2.39 (1.5cr) and -2.56 < [(0/H)mw] < -2.10 {2a). 
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We conclude that the LCP measurements constrain the 
mass-weighted mean oxygen abundance of the IGM to 
be [(0/H)mw] = -2.45 ±0.06 (I.Sct) and [(0/H)mw] = 
—2.33 ± 0.23 {2a). We can compare this result directly 
with the mass- weighted mean in AOS of [(0/H)mw] = 
— 2.14±0.14 (Icr), which is marginally consistent with our 
results. The volume- weighted mean oxygen abundance is 
more weakly constrained to [(0/H)vw] — —3.06 ± 0.22 
(1.5cr) and [(0/H)vw] -2.74 ± 0.64 (2cr). These 
error estimates do not include systematic uncertainties 
associated with the shape and intensity of the ioniz- 
ing background (see §6.ip or uncertainties in th e so- 
lar oxygen abundance. Using ICaffau et al.l (|2008l ) val- 
ues for the solar abundance of oxygen would give a 
mass-weighted [(0/H)mw] = -2.34 ± 0.06 (1.5cr) and 
[(0/H)mw] = -2.22 ± 0.23 (2a). 

With this measure of [(0/H)mw]) and extrapolating to 
all densities, we can calculate the cosmic oxygen density 
using 

f^Oxy,iGM « Y^Aon,l0i(^/^)Mw] (0/H)q, (7) 

where Yh — 0.76 is the mass fraction of hydrogen and 
Aq = 16 is the atomic mass of oxygen. We obtain 
a value of floxy.iGM = 1.4(±0.2) x lO^^ (l.5cr). This 
would be ~ 5% smaller if we limited ourselves to the 
density range — 0.5 < \og[p/p] < 2 (the range for which 
the mass- weighted mean was calculated). This shows 
that the calculation is not sensitive to the details of the 
density range integrated over, but this is not a state- 
ment about the metal enrichment of systems beyond our 
constraining density range. In particular, the popula- 
tion of strong absorbers, such as those resulting in the 
Thy a ~ 3.3 point of our LCP, are b etter dealt with in 
directed searches (iFrank et"alll2009a l |bh . 

If we follow the procedure of iBouche et al.l (|2007( ) (and 
AOS) and take Ozjgm ~ ^^Oxy,iGM/0.6, then we ob- 
tain a total cosmic metal density of the IGM, J^zjom ~ 
2.3 X 10-6. This is w 16% of their metal budget of 
flz — 1.5 X 10^^ at z = 3 (from their figure 1 and based on 
integrated cosmic star formation with standard assump- 
tions about about stellar IMF and yields). This metal 
content is over half that found in galaxies at z = 2.5 
(30%) wh ile sub-DLAs m ay contribute between 2 and 
17% (Bouche et al.ll2007[ ). Summing these values would 
account for over half of the metals produced by star for- 
mation. The remainder could plausibly resid e in interme- 
diate d ensities (e.g. strong O VI absorbers, iFrank et al.l 
l2009bD or in a hidden warm-hot phase, but uncertain- 
ties in the calculated metal production and the observa- 
tional estimates preclude strong statements about "miss- 
ing" metals. 

6. DISCUSSION 

6.1. Ionization Corrections 

We have obtained the most precise constraints to date 
on the mass- weighted mean O/H and iloxy, but our anal- 
ysis rests on the assumption that the detected O VI is 
mostly ionized by a UV background with the spectral 
shape of the HMOl QG model. The use of this UV back- 
ground model is motivated by the work of S03 and AOS. 
They also use two other background models: one result- 
ing from quasars only ("Q") and a softened version of 
the QG model ("QGS") with a lower flux above 4 Ryd. 



TABLE 1 

Mass- WEIGHTED oxygen abundance in some best fitting 

MODELS 



Model (fO/Hl, (p/p)cut) 


[{0/H)mw1 


Goodness-of-fit 


-2.15, 


-2.15 


1.8(T 


-2.10, 1 


-2.25 


1.8cr 


-2.00, 2 


-2.39 


1.5cr 


-1.80, 3 


-2.41 


1.3cr 


-1.65, 4 


-2.45 


1.2cr 


-1.45, 5 


-2.43 


1.3cr 


-1.80, 6 


-2.43 


1.4cr 


-1.15, 7 


-2.43 


1.7a 


-0.95, 8. 


-2.36 


l.Sa 


AOS 


-2.36 


2a 



Note. — For models along the ridge- line of Figure [T] with en- 
richment density cutoffs but no scatter in metallicity at fixed den- 
sity, we list the [O/H] corresponding to the mass- weighted mean 
oxygen abundance ratio (O/H) in column 2 and the acceptability 
of the fit to the data in column 3. 

They find that those two models are inconsistent with 
their relative abundance measures of [0/C] and [0/Si], 
while they obtain a good fit using the QG model. It 
is not clear to what extent small variations in the QG 
model can be tolerated, and what impact these would 
have on the quoted results. This is clearly an area that 
requires further investigation. 

Our simulations assume a power-law equation of state 
and so under-produce gas at T > lO^K, which is domi- 
nated by coUisional ionizatio n. This cho i ce is justified in 
light o f stu dies by|C arswell et al.l (|2002[ ). iBergeron et al.l 
(|2002t l and lSimcoe et al. f)200l who show that O VI lines 
are typically too narrow to be in this phase. It is also no- 
table that O VI which arises from collisional ionization 
also tends to be strong due to the high O VI fraction 
in this regime, while the O VI absorption probed here 
is relatively weak (even compared to the studies listed 
above) . 

AOS find a substantial increase in [O/H], for a fixed 
[0/C] ratio, for systems that have p/p > 30. Lower den- 
sities are well described by photoionization and a fixed 
[0/C] ratio in their models. They conclude that p/p > 30 
regions are in fact coUisionally io nized (while systems 
with p/p < 30 are photoionized). ISimcoe et al.l (|2004[ ) 
also find such a change in [0/C] and infer a change in 
ionization mechanism at higher densities. We find indi- 
cations of the same effect for our measure at T^ya ~ 3.3, 
which shows O VI absorption at a level far higher than 
for lower TLya- We have omitted this data point for our 
photoionization dominated analysis. We have only one 
point showing this effect, since we do not use higher or- 
der Lyman lines in our analysis and are therefore unable 
to cleanly isolate the highest density regions. We find 
a good fit at all detectable Lya absorption levels below 
Thya ~ 3.3 using a constant [O/H] and assuming a pho- 
toionized medium. 

One should note that AOS normalise their UV back- 
ground to a level that differs from ours. This leads to 
an adjustment on the ionization correction. The ratio of 
optical depths is related to the metallicity by 



" T OVI ■ 


= log 


"OVI 




L Thya J 




■ no 


nm- 



log^ ^log™l!£l +[0/H]+log(0/H)o + 5, 

(8) 

inhere S = log[(/A)ovi/(/A)LyQ]- The first term on the 
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Fig. 9. — The ionization correction parameter 

("OVI'^h)/('^0'^hi) in terms of the UV background inten- 
sity (normalized to the HMOl level) and ionization rate at 
z = 3. This is calculated for 3 densities that span our range of 
sensitivity: mean density {solid line), log[p/p] = 0.5 {dotted line), 
and log[p/p] = 1 {dashed line). The thin vertical lines show the 
UVB intensity used in AOS and here. 

right hand side of this equation is the ionization correc- 
tion term, and this is plotted against UV background 
intensity (normalized to the HMOl level) and the ion- 
ization rate in Figure [HI This was calculated assuming 
the equation of state in equation |3] and z = 3 for three 
overdensities: \og[p/p] = 0, 0.5 and 1 that span our sen- 
sitivity range. Our chosen UV background intensity of 
1.35 X HMOl (r = 13 X 10-"s-i) is shown along with 
the AOS value of 0.48 x HMOl (F = 4.6 x lO-i^g-i). 
For mean density the ionization correction is largely un- 
changed between our simulations and those of AOS, how- 
ever, our models would result in systematically lower 
[0/H] resulting from oxygen in systems of log[p/p] > 0.5 
by between 0.3 and 1.1 dex. This might go some way to 
explaining the 0.2 dex lower value we find for [0/H]. An- 
other source of potential discrepancy is differences in the 
density distributions between the studies. Despite these 
issues, we are in broad agreement with a metal budget 
fraction of - 15 - 35% in AOS for z = 2 - 3. 

The success of the LCP search method for detection of 
O VI in the SDSS spectra may indicate that it is appeal- 
ing for other searches for low level absorption in the IGM, 
where large numbers of moderate resolution spectra with 
moderate signal-to-noise and a substantial uncertainty 
in the continuum are available. There are a number of 
other species that may be detectable in the SDSS sample 
with a similar approach. Correlation searches in spectra 
from other instruments such as FUSE and STIS, or even 
HIRES and UVES, may benefit from this kind of analy- 
sis, which is resistant to errors in the continuum. Mod- 
ulation of the number of local pixels used may also be a 
useful approach to measuring clustering of absorbers. 

6.2. Continuum Correction Uncertainty 

We systematically underestimate our continuum fit- 
ting such that the continua must be corrected up (and 
the transmitted flux down) by 16% to match the mean 
flux measured in studies using high resolution spectra. 
The uncertainty in those measurements may lead to sys- 



tematic errors in our me asurements. We use (F ) = 
0.68 ± 0.04 atz = 3 from iMcDonald et al.l (|2006l ) and 
those la error estimates correspond to allowed correc- 
tions of 11-25%. This would lead to a 15% offset in 
{T-Lya), which would be corrected by a 15% offset in the 
recombination rate. This change in recombination rate 
results in change to the H I fraction which is indepen- 
dent of density for the regime probed here, leaving the 
Lya forest largely unchanged. There is a small distortion 
due to the fact that an error which is a multiple in flux 
would be corrected by a factor which is a multiple in op- 
tical depth, but crucially most of the signal in the LCF 
search comes from higher optical depths for which the 
correction is even smaller. Also large errors of this sort 
would be evident in the comparison between the FDF of 
the flux in the simulations and the data in Figure IH 

A systematic error in the chosen F (based on an error 
in the mean flux in the literature) would also impact 
upon the O VI fraction. This error corresponds to half a 
minor tick interval in the upper horizontal axis of Figure 
[5] and so ionization correction parameter would be largely 
unchanged over the density range of interest. One can 
conclude that in the context of (related) uncertainties in 
the shape and amplitude of the UV background and so F, 
these issues are a sub-dominant systematic uncertainty. 

6.3. Other Searches for O VI in the ICM 

iFox et all (|2008h argue that the lack of velocity de- 
pendence in the column densities {N) of their proximate 
{5v < 8000 kms~^ from zqso) sample indicates that 
their 26 weak O VI absorbers must be coUisionally ion- 
ized. The majority of these absorbers likely correspond 
to systems of p/p > 10 based on their associated H I col- 
umn d ensities and using the A^— p/p relation from lSchayd 
(|2001| ). Hence they are sampling predominantly stronger 
absorbers than those considered here, so this work is not 
at odds with our choice of a photoionization model for 
our calculation of [0/H]. 

Our results appear consistent with other comparable 
measures of [0/H] in the IGM, given the differen ces of 
metho d and thus the systems sampled. Cars wel l et alJ 
(|2002f ) also use high resolu t ion sp ectra to search for O VI, 
but unlike lAguirre et al.l ()2008[ ) they use Voigt profile 
fitting of detected lines. They argue that these lines at 
z = 2 are photoionized (and so are comparable with those 
w e find), wit h abundances 10""^ — 10~^ of solar. 

iTelfer et al.l (|2002l ) use stacked O V in FUSE spectra 
to determine [0/H] and obtain a range of -2.2 to -1.3. 
Their higher value is partly caused by a combination of 
their use of a "Q" UV background model and by their 
lower redshift range (1.6 < z < 2.9). The main cause for 
the discrepancy is likely their use of the 78th percentile of 
their stack (in ord er to maximiz e signa l-to-noise), while 
we use the 50th. lAguirre et all ()2008f ) argue that this 
difference likely correspon ds to a 0.5 dex adjustment, 
bringing the results of Telfer et al.l (2002) broadly into 
li ne with ours. 

Her geron fc Herbert-FortI (|2005f l search for O VI lines 
at somewhat lower redshift range (2 < z < 2.6) than 
ours. They divide these lines as 'metal-poor' and 'metal- 
rich' . This differs greatly from our analysis and so is 
difficult to compare. In particular they find the [0/H] 
of systems with detected O VI absorption while we cal- 
culate the average metallicity over all Lya absorption. 
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• A model with constant [0/H] = -2.15to:o9 P'^O" 
vides an adequate fit to the data (at the 2a level), 
but a fit can also be obtained by removing metals 
from lower density regions and elevating the metal- 
licity in regions enriched. 

• The addition of lognormal metallicity scatter leaves 
the quality of fit of these models largely unchanged, 
where the scatter free model is treated as a measure 
of the mean 0/H (or (Z)) and not the mean [0/H] 
(or (logZ)). 

• The mass-weighted mean oxygen abundance is 
nearly constant among our viable models. As 
a result we have been able to place the tight- 
est constraints thus far on this quantity, with 
[(0/H)mw] = -2.45 ± 0.06 (at 1.5a confidence) 
computed over the density range —0 .5 < l og[p/p] < 
2. This value is based on a iGrevessc ct al] 



1996 ) solar abundance of oxygen. A ICaffau et al 



20081 ) solar abundance would give a abundance 



Their 'metal-poor' sample shows —3 < [O/H] < —1 and 1 p/p ^ 9 s-nd can conclude the following: 
a cosmic oxygen density fioxy ~ 2.3(±Q;g) x 10~^ using 
a photoionization model similar to the "Q" model de- 
scribed above, which is consistent with our results, but 
t he numbers a r e not d irectly comparable. 

iFrank et all ()2009a| ) have searched the same SDSS 
dataset with the goal of identifying strong O VI ab- 
sorbers, requiring additionally associated Lya and Ly/? 
absorption. The approach used is complementary to this 
work, as they can only retrieve strong absorbers that are 
typically saturated in both these Lyman lines. Hence, 
these absorbers are a different population from the ones 
probed her e with higher density and (perhaps) collisional 
ionization. iFrank et al.l (|2009bf ) provide a measure of the 
cosmic oxygen density contributed by these systems of 
f^ovi > 1.9 X 10~^h~^ (which corresponds to a firm lower 
limit of rioxy > 1 X assuming the maximum 

O VI ionization fraction of 20%) at absorber redshifts 
2.8 < Zabs < 3.2. 

Our measure of the fraction of the estimated metal 
bu dget detected at z = 3 is in line with values calculated 
bv lBouche et al.l (|2007t ) using results from various papers 
and various metals. They infer a 'forest' contribution of 
< 15% at z = 2. One study they quote is particularly no- 
table: ISimcoe et al] (|2004[ ) find fl pxy » 2 x 10"^ using 
a 'hard' ('Q' like) UV background. iBouche et all ()20071 ) 
reanalyze these results in the context of a 'softer' ('QG' 
like) UV background model and find fioxy ~ 5.0 x 10~^ 
at z = 2.5, and so a fraction of the metal budget of 
« 30%. This measure is consistent with our findings 
given redshift evolution and systematic uncertainties as- 
sociated with search methods and ionization corrections. 
ISimcoe et aP (|2004l ) find a mean [0/H] = -2.85 for their 
lines (which is effectively a volume-weighted mean) for 
the 70% of systems in which they detected oxygen, which 
is also consistent with our findings (setting the other 30% 
to zero metallicity would result in a volume-weighted 
mean of [(0/H)vw] = -3.00). 

7. CONCLUSIONS 

We have developed a new method to detect O VI ab- 
sorption in low-density regions with a data set that might 
initially seem ill suited to this purpose, the moderate res- 
olution {R — 1800), moderate signal-to- noise ratio SDSS 
spectra of high-redshift QSOs. Not only have we success- 
fully detected O VI with high significance (the null case 
is ruled out with reduced = 10), but we have placed 
tighter limits on the mass-weighted mean oxygen abun- 
dance of low-density regions than ever before. The suc- 
cess of this method hinges on the sheer size of the SDSS 
sample of spectra and the use of local, uncorrelated pixels 
near the O VI absorber as a way of characterizing both 
the continuum fitting error and the degree of contami- 
nating absorption. This locally calibrated pixel (LCP) 
search for O VI seen in SDSS DR3 is the main result of 
the paper and the high precision of this new tracer for 
metals in the IGM is clear. Modeling has been necessary 
to interpret these results and conclusions about oxygen 
abundance the lognormal random fields have been used 
to produce more than 4.5 million synthetic spectra. 

Our search at 2.7 < z < 3.2 is not sensitive to the 
density dependence of the [0/H] because of spectral res- 
olution constraints, but using detailed modeling we are 
able to measure the oxygen abundance IGM densities of 



[(0/H)mw] = -2.34 ±0.06. 

• The models that fit the data have volume- weighted 
mean, [(0/H)vw] = -3.01 ± 0.33 (at 1.5a confi- 
dence) . 

• We calculate the cosmic density of oxygen, 
f^OxyjGM = 1.4(±0.2) X l O"*". In the context of 
the work bv lBouche et al.l (|2007t ) this constitutes a 
total metal contribution of fizjcM ~ 2.3 x 10~^, 
which is « 16% of their estimated metal budget at 
at z = 3. 

This novel LCP search method need not be limited 
to searches for metals in the IGM. It could be adapted 
to spectral searches wherever bulk, low-level correlated 
absorption is expected and continuum errors and con- 
taminating absorbers are the limiting factor. Little prior 
knowledge of these errors and uncertainties is required, 
only that they are uncorrelated. 
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